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Abstract 

Using an effective interaction approach to describe the interactions between the dark matter 
particle and the light degrees of freedom of the standard model, we calculate the gamma-ray flux 
due to the annihilation of the dark matter into quarks, followed by fragmentation into neutral pions 
which subsequently decay into photons. By comparison to the mid-latitude data released from the 
Fermi-LAT experiment, we obtain useful constraints on the size of the effective interactions and 
they are found to be comparable to those deduced from collider, gamma-ray line and anti-matter 
search experiments. However, the two operators induced by scalar and vector exchange among 
fermionic dark matter and light quarks that contribute to spin-independent cross sections are 
constrained more stringently by the recent XENON100 data. 



1 



I. INTRODUCTION 



The presence of cold dark matter (CDM) in our Universe is now well established by a 
number of observational experiments, especially the very precise measurement of the cosmic 
microwave background radiation in the Wilkinson Microwave Anisotropy Probe (WMAP) 
experiment The measured value of the CDM relic density is 

n CDM h 2 = 0.1126 ±0.0036 , (1) 

where h is the Hubble constant in units of 100 km/Mpc/s. Though the gravitation nature 
of the dark matter is commonly believed to be well established, its particle nature remains 
allure except that it is nonbaryonic and to a high extent electrically neutral. 

One of the most appealing and natural CDM particle candidates is the weakly-interacting 
massive particle (WIMP). It is a coincidence that if the dark matter (DM) \ is thermally 
produced in the early Universe, the required annihilation cross section is right at the order 
of weak interaction. The relation between the fractional relic density of x relative to the 
critical density and its thermal annihilation cross section can be given by the following simple 
formula ^ 

Olpb 

U x h _ ^ , [2) 

with (av ) being the annihilation cross section of the dark matter around the time of freeze- 
out, at which the annihilation rate could no longer catch up with the Hubble expansion rate 
of the Universe. Assuming the measured QcBMh 2 to be saturated by a single component 
WIMP, its annihilation cross section should be about 1 pb or 3 x 10 -26 cm 3 s _1 . This 
is exactly the size of the cross section that one expects from a weak interaction process, 
which implies an appreciable size of production rate of the WIMP at the Large Hadron 
Collider (LHC) as well as the event rates for direct and indirect searches that reach the 
sensitivities of dark matter experiments like XENON100 and Fermi-LAT respectively. In 
general, production of dark matter at the LHC would give rise to a large missing energy. 
Thus, the anticipated signature in the final state is high-p-r jets or leptons plus a large 
missing energy. Note that there could be non-thermal sources for the dark matter, such as 
decay from exotic relics like moduli fields, cosmic strings, etc. In such cases, the annihilation 
cross section in Eq. (T5]) can be larger than the value quoted above. 

There have been many proposed candidates for the dark matter. Without committing 



to any particular DM model so as to perform a model independent analysis, we adopt an 
effective interaction approach to describe the interactions of the dark matter particle with 
the standard model (SM) particles J3 - 15] - One simple realization of the effective interaction 
approach is that the dark matter particle exists in a hidden sector, which communicates to 
the SM sector via a heavy degree of freedom in the connector sector. At energy scale well 
below this heavy mediator the interactions can be conveniently described by a set of effective 
interactions. The strength of each interaction depends on the nature of the dark matter 
particle and the mediator. The most important set of interactions are among the fermionic 
dark matter \ and the light quarks q described by the effective operators (xTx) iQ^'q) where 
T and r" are general Dirac matrices contracted with appropriate Lorentz indices. We will 
discuss these and other operators in more details in the next section. 

There have been some recent works on constraining the interactions at present and future 
collider experiments J* 



search experiments [13 



i |6|, using gamma-ray experiments [7 10 , [ifil, and using anti-matter 
15[. There was another work in which the dark matter couples only 
to the top quark and corresponding predictions at direct and indirect detection experiments 



as well as colliders were obtained 



12] . It was also shown in Ref . 



that additional radiation 



of electroweak bosons in the final state can modify the energy spectrum, especially at the 
lower end of the spectrum. Lifting the helicity suppression by radiating off an electroweak 
gauge boson from the external light fermion legs due to Majorana dark matter annihilation 
is emphasized in Ref. [ji], 19]. 

In Ref. js], monochromatic photon-line flux was calculated via a loop with quarks running 
in it and photons being attached to the internal quark line. Although the photon-line would 
be a smoking-gun signal to compare with the data, the rate is suppressed because of the 
loop factor. On the other hand, photons can come from the decay of neutral pions, which 
in turn come from the fragmentation of the quarks in the annihilation of the dark matter. 
The chance that an energetic quark fragments into neutral pions is high and the branching 
ratio of a neutral pion into two photons is 98.823% 20|. Therefore, the amount of photons 
coming from the quark fragmentation is much larger than those coming off a loop process. 
Nevertheless, the spectrum of such photons is continuous and in general have no structure, 
except for a cutoff due to the mass of the dark matter. In this work, we focus on the 
continuous gamma-ray flux spectrum coming from the fragmentation of quarks into neutral 
pions, followed by their decays into photons, in the annihilation of the dark matter. Such 



annihilation of the dark matter will give rise to an additional source of diffuse gamma-rays 
other than the known backgrounds. If the experimental measurement is consistent with 
the known gamma-ray background estimation, then one could use the data to constrain the 
amount of gamma-ray flux coming from the dark matter annihilation, thus constraining the 
effective interactions between the dark matter and the quarks. 

The data on the photon spectrum from the mid-latitude (10° < \b\ < 20°, 0° < I < 360°) 



2l| recorded by the Fermi-LAT indicated a continuous spectrum and mostly consistent 
with the known backgrounds. We can therefore use the data to constrain on additional 
sources of gamma-ray, namely, the annihilation of the dark matter into quarks, followed 
by fragmentation into neutral pions, which further decay into photons. The production 
of photons via neutral pions is the dominant mechanism for gamma-rays. There are also 
other ways that the quarks from annihilation of dark matter can produce gamma rays, such 
as bremsstrahlung off quark legs, synchrotron radiation, or inverse Compton scattering on 
background photons, but these processes are all a cm suppressed relative to the fragmentation 
of the quarks into neutral pions. We focus on the fragmentation of the quarks coming from 
the annihilation of the dark matter as the signal in our analysis. We employ two approaches 
of obtaining the photon spectrum due to fragmentation of light quarks, (i) We use the 



process e + e~ — > qq with initial radiations turned off in Pythia 22] and extract the photon 
spectrum in the final state. The photon mainly comes from the decay of 7r°, which are 
in turn produced by fragmentation of light quarks, plus a very small fraction from the 
bremsstrahlung photon off the quar k legs, (ii) We use the fragmentation function of q, q, g 



into 7T° from the fitting of Ref. 23|, then convolute with the dN / dE^ir — > 77) to obtain 
the photon spectrum of quarks into photon. We found that both approaches give almost the 
same photon spectrum from light quarks. The resulting limits using both approaches are 
also the same within numerical accuracy. However, we have to use the second approach when 
we place limits on the operators involving gluons, because we do not find an appropriate 
process in Pythia for extraction of g — > 7 fragmentation. 

The choice of the mid-latitude data instead of the Galactic Center is simply because the 
gamma-ray in this region is dominated by local sources and we have clarity in understanding 
the background flux and point sources within the mid-latitude. On the other hand, the 
Galactic center is supposed to have a number of known and known-unknown point sources, 
including a supermassive black hole near the Center, and perhaps some unknown sources too. 



Given the purpose of constraining the new DM interactions it is better to understand clearly 
about the background in mid-latitude, rather than the larger flux from the Galactic Center. 
The Galactic diffuse gamma rays originate primarily from the interactions of high energy 
charged particles contained in cosmic rays with the nuclei in the interstellar medium and the 
associated radiation fields of the charged particles, via a few mechanisms briefly described 
in Sec. III. While most of them are well understood, the extra-galactic component has a 
larger uncertainty. We will choose a normalization such that the total background diffuse 
gamma-ray flux is consistent with the Fermi-LAT measurement of diffuse gamma-ray flux 
in the mid-latitude. This approach is the same as the Fermi-LAT when they estimated the 



extra-galactic diffuse component [21 ]. 



The organization of the paper is as follows. In the next section, we describe the inter- 
actions between the dark matter particle and the SM particles, in particular quarks and 
gluons. In Sec. Ill, we discuss various sources of diffuse gamma-ray flux that constitute the 
known background and calculate the gamma-ray flux due to the dark matter annihilation 
using the effective interactions. We compare with other constraints and conclude in Sec. IV. 

II. EFFECTIVE INTERACTIONS 

For simplicity, we will assume there is only one component of dark matter denoted by x 
and it is a standard model singlet. 

The first set of operators we will be considering is for fermionic DM and its effective 
interactions with light quarks via a (axial) vector- or tensor-type exchange are given by the 
following dimension 6 operators 

A=i- 6 = Sr (WX) (qT'q) , (3) 

i 

where T, V = 7 M , 7^7 5 , a^ u or a^ v rf with a^ iu = 2(7^7^ — 7 1/ 7 /i )/2. A* is the heavy scale for 
the connector sector that has been integrated out and C is an effective coupling constant of 
order 0(1). It is understood that for Majorana fermion the vector and tensor structures of V 
are absent. Thus, for vector or tensor type interaction the fermion \ in Eq.([3]) is understood 
to be Dirac. We will be focusing on Dirac fermionic DM in this work, but our results are 
also applicable to Majorana dark matter. Note also that due to the following identity 

c^ 7 5 = \e^a aP , (4) 



TABLE I. The list of effective interactions between the dark matter and the light degrees of 



freedom (quark or gluon). We have suppressed the co 



These operators have also been analyzed in Refs 



. s,y,y,Q]. 



or index on the quark and gluon fields. 



Operator Coefficient 


Velocity Scaling in (av) 


Dirac DM, (axial) vector/tensor exchange 


0\ = (xi^x) (qi^q) 

2 = (ry'Vx) (qi^q) 

03 = (x7 M x) (qi^q) 

04 = (x7 M 7 5 x) (?7ai7 5 <?) 

05 = (xa^x) (q^uq) 
G = (x^Vx) (q^q) 


c 

A? 

c 
c 
a 

G 

c 

A 7 


2 2 
mtv 

m 2 x 

2 2 
m^v 

m 2 x 
m x 


Dirac DM, (pseudo) scalar exchange 


o 7 = (xx) (qq) 
Os = (x7 5 x) (qq) 
O9 = (xx) (qj 5 q) 
do = (x7 5 x) (qi 5 q) 


Cm q 

A» 
iCrriq 

A^ 
iCrriq 

A 3 
Cm q 

A a 


222 
m^m^v 

2 2 
m q m^ 

222 
m*mtv 

2 2 
m q m^ 


Dirac DM, gluonic 


On = (XX) G^CT 

012 = (x7 5 x) G^G^ 

013 = (XX) G^G^ 
Ou = (xi 5 x)G^ 


Ca a 
4A 3 " 

iCa s 
4A a 

Ca a 
4A 3 " 

iCa s 
4A 3 


4 2 
m^y 

4 

m x 

4 2 
mtv 

A 

4 

m x 


Complex Scalar DM, (axial) vector exchange 


Oi 5 = (x ] d^x) {rfq) 
Oie = (x f <Cx) (97"7 5 9) 


C 
A 7 

C 
A 7 


2 2 

A 

2 2 

mtv 

A 


Complex Scalar DM, (pseudo) scalar exchan^ 


;e 




On = (x ] x) (qq) 
Oie = (xh) (qi b q) 


Crriq 

iCm q 
A y 




Complex Scalar DM, gluonic 


019 = (x f x) G^ 

20 = (x t x)0^0^ 


C« s 
4A 7 

4A 2 


m x 
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the axial tensor cr' JiUr y 5 is related to the tensor a a ^ and thus should not be regarded as an 
independent set for Dirac fermionic DM. However, for Majorana fermionic Xi such axial 
tensor structure can be present. 

Next set of operators are associated with (pseudo) scalar-type exchange 

4=7-io = ^ (xT X ) (qT'q) , (5) 

i 

where r, V — 1 or ry 5 . The m q dependence in the coupling strength is included for scalar- 
type interactions in accord with the trace anomaly in QCD. We use the current quark masses 



in the Lagrangian given by 



20j: 



m u = 0.0025 GeV, m d = 0.005 GeV, m s = 0.101 GeV, 
m c = 1.27 GeV, m b = 4.19 GeV, m t = 172 GeV. 

Another light degree of freedom that couples to the Dirac dark matter is the gluon field 







£,.„-i2 = CC "}* nd (XT X ) G""Gl„ (6) 

i 

4=is- M = ^f.?^ (XT X ) G^Gl„ (7) 

i 

where T = 1 or 27 s . For operators involving gluons, the factor of strong coupling constant 
a s (2m x ) is also included in accord with the trace anomaly and is evaluated at the scale 2m x 
where m x is the dark matter mass. 

Finally, we also write down the corresponding operators for complex scalar dark matter. 
Again, we note that the interactions for real scalar dark matter is similar to complex one 
and differ by a factor of two. We simply focus on the complex scalar dark matter. The 
operators corresponding to vector boson exchange are 

4=16,16 = § (x^Xx) (g7 M r g ) , (8) 

where T = 1 or 7 5 and x^Xx — xK^^x) ~ (dfiX^)x- Those corresponding to a scalar boson 
exchange are 

4=i 7 ,i8 = ^ u f x) m) , (9) 



1 We do not study the other gauge bosons, like W and Z bosons, because they decay into light quarks 
which then fragment into photons, would be softer in this case. 
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where T = 1 or 27 s . The corresponding gluonic operators are 

4=19 = Ca fj? x) (xh) G^G a , v , (10) 

i 

C i=w = lCa ^ ( x t % ) G ^G ailv . (11) 
j 

The whole set of operators we are studying are tabulated in Tabled] together with their corre- 
sponding coefficients and scaling with velocities in the annihilation cross sections. Without 
a particular model in mind we will treat each interaction independently in our analysis by 
considering one operator at a time and setting the coefficient C = 1 for simplicity. □ The 
relative importance of each operator can be understood by considering the non-relativistic 
expansion of the operator and studying the velocity dependence. It was fully discussed in 
Ref. 13] and hence we only briefly summarize here for convenience. In the non-relativistic 
limit, the spinors for the Dirac DM x an d X annihilation are ip ~ (£, e£) T and z/i ~ (er^, 7^)7° 
where £ and 77 are two-components Pauli spinors and e = 0(v/c). We can expand %jj^^ip as 

^7> ~ (1 + e 2 W°ii 

where the spatial components are not suppressed by v/c. On the other hand, ip^^^ip in the 
non-relativistic limit are 

~ (1 + e 2 )^ 

where the spatial components are now suppressed by v/c. It is clear that in the non- 
relativistic limit the time and spatial components of the vector and axial vector bilinear 
behave very differently. We can then consider them separately when it is contracted with the 
trace of the light quark leg. If we look at the trace of (qj^q) or (qj^^q) in the annihilation 
amplitude, the time component part after being squared gives a quantity close to zero, while 
the spatial component part gives a quantity in the order of m 2 . Therefore, it is clear now 
that ip-f^ip multiplied to (qj^q) or (qj^q) will not be suppressed, while t/j^^ip multiplied 
to (qj/j.q) or (qj^q) will always be suppressed. Therefore, the operators 0\ and O3 can 
contribute to annihilation much more than the operators O2 and O4. Thus, the limits on 



For gluonic operators, their coefficients Cs are induced at loop level and can be smaller. 
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0\ and O3 are much stronger than O2 and O4. All the other operators listed in Table H] can 



be understood similarly 



131 ] . We note that some of the operators are doubly suppressed by 



the velocity of the dark matter combined with either a light quark mass or strong coupling 
constant. 

In the calculation of the gamma-ray flux from dark matter annihilation presented in the 
next section, we only include the light-quark flavors. We ignore the XX ~^ contribu- 
tion, because the t and t first decay into bW —> bqq' before each light quark undergoes 
fragmentation into hadrons, including pions. Therefore, the gamma-ray spectrum would be 
significantly softer than the direct fragmentation as in XX ~* QQ 13]. 

Note that some of the lower limits that we obtain in Table [TTT1 are relative low compared to 
the dark matter mass. Notably for the operators O7 9, On t i 3 , Oi 5j i 6 and Oi 7j i 8 . In such cases, 
one may question the validity of the effective interaction approach. The physics behind is 
easy to understand. The effects of such operators are very suppressed because of the small 
velocity suppression or helicity suppression, not because of the size of the A. Therefore, the 
A has to be small enough in order to see an effect from these operators. We argue that 
the effective momentum transfer of such velocity-suppressed operators should be m x (v/c). 
With (v/c) ~ 10~ 3 for the DM velocity at the present epoch, as long as the ratio m x (v/c)/A 
remains small, we expect the effective interaction approach can still be valid. 



III. GAMMA-RAY FLUX 

A. Background Diffuse Gamma Rays 

The Galactic diffuse gamma rays originate primarily from the interactions of high energy 
charged particles contained in cosmic rays with the nuclei in the interstellar medium and the 
associated radiation fields of the charged particles, via a few of the following mechanisms. 

(i) Gamma-rays coming from the ir° decay, which in turn comes from the interactions 
of the cosmic rays with the nucleons in the interstellar medium. This dominates the 
background flux for energy higher than 1 GeV. 

(ii) Inverse Compton scattering occurs when high energy e ± collide with the photons of 
the interstellar medium. 
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(iii) Bremsstrahlung photons occur when high energy are deflected by the Coulomb field 
of the interstellar medium. 

(iv) Synchrotron radiation occurs when high energy e ± are deflected by Galactic magnetic 
field. 

(v) An extragalactic background, which is expected to be isotropic and receives contri- 
butions from many sources including unresolved point sources, diffuse emission from 
large scale structure formation and from interactions between ultra-high energy cosmic 
rays and relic photons, etc. This background is the least determined and so a fairly 
large uncertainty is associated with it. The Fermi-LAT has a measurement of diffuse 
gamma-ray in the mid-latitude region and fitted the extra-galactic background by 

where A and 5 are fitted parameters (the power- law index is 7 = 1 8 — 2 1 ) . In Ref . 21] , 
the power-law is fitted to be 7 = 2.41 ± 0.05 and A can be determined by the to- 
tal flux of EGB ("extragalactic" diffuse gamma-ray emission) as A = (0.95 io'17) x 
10~ 6 GeV cm~ 2 s" 1 sr" 1 for E > 100 MeV. 

Since Fermi-LAT did not give details about the parameters of GALPROP (Cosmic ray 



propagation code) |24j that they used in their treatment, we run GALPROP (the web 
version) to obtain the various diffuse Galactic backgrounds (i) to (iv) and fit EGB component 
in Eq. ffl2l to Fermi-LAT data. The relevant GALPROP parameters that we used are shown 
in Table [Til Our fitted EGB is given by 

r!(h / VP \ ~0.36 

E2 IE = <°" * 10 ~ 6) {lusTGev) GeV Cm " S " (13 > 
which gives a power-law index 7 = 2.36 and E is in GeV. It is close enough to the one 
obtained by Fermi-LAT. The x 2 = 2.435 for 7 d.o.f. The various curves are all within 
the uncertainties quoted in the Fig. 6(a) of the Fermi-LAT paper 21|. Various diffuse 
background curves and their sum are shown in Fig. [TJ 



B. Dark Matter Annihilation 



The dominant DM contribution to photon flux in this scenario comes from 

XX ->• qq ->• tt° + X -> 2 7 + X , 
10 



(14) 




10 -7 I , ■ ■ x ■ ■ ■ 

1 10 100 

Ey(GeV) 

FIG. 1. The photon spectrum E 2 (d^/dE 1 ) versus the photon energy for the diffuse background 
gamma rays, including ir° decays (short-dashed-dotted) from interactions of cosmic rays with in- 
terstellar medium, inverse Compton scattering IC (long-dashed-dotted) , bremsstrahlung (dotted) 
and extragalactic EGB (short-dashed). Their sum is shown as the lower solid line. A dark matter 
component due to DM annihilation of XX ~~ QQ is a l so shown (long-dashed) and added to the total 
background (upper solid line). The DM annihilation is due to a 200 GeV DM particle with the 
effective interaction operator 0\ and A = 1.5 TeV. 



in which all the q,q (q — u, d, c, s, b) have probabilities fragmenting into n°, which then 
decay almost entirely into two photons. As mentioned in the Introduction, we employ two 
approaches of obtaining the photon spectrum due to fragmentation of light quarks, (i) Using 



the process e + e~ — > qq with initial radiations turned off in Pythia [22|] and extracting the 
photon spectrum in the final state. The photon mainly comes from the decay of tt°, which 
are in turn produced by fragmentation of light quarks, plus a very small fraction from the 
bremsstrahlung photon off the quark legs, (ii) Using the fragmentation function of q, q, g 
into 7T° from the fitting of Ref. 23j, then convoluting with the dN / dE 1 (ir° — > 77) to obtain 
the photon spectrum of quarks into photon. Both approaches give the photon spectra close 
enough to each other for our purpose of numerical calculations. Thus, the resulting limits 
using both approaches are also the same within numerical accuracy. However, we have 
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10 -7 I , , ^ J , : \ 

1 10 100 

Ey(GeV) 

FIG. 2. Same as Fig. [I] but the DM annihilation is due to a 50 GeV dark matter particle with the 
effective interaction operator 0\ and A = 0.87 TeV. 



10~ 5 




10"' 



10 



EJGeV) 



100 



FIG. 3. Same as Fig. Q] but the DM annihilation is due to a 500 GeV dark matter particle with 
the effective interaction operator 0\ and A = 1.9 TeV. 



to use the second approach when we place limits on the operators involving gluons in the 
annihilation, because we do not find a simple process in Pythia for extraction of g — > 7 
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TABLE II. The values of the parameters that we used in running the code GALPROP to reproduce 
the background curves in Figs. [TH3J For unstated parameters we employ the default values. 



Parameter (unit) Value 

Minimum galactocentric radius r m j n (kpc) 00.0 

Maximum galactocentric radius r max (kpc) 25.0 

Minimum height z m { n (kpc) -04.0 

Maximum height z max (kpc) +04.0 
ISRF factors for IC calculation: optical, FIR, CMB ISRF factors 1.9, 1.9, 1.9 

Gamma-ray Intensity Skymap Longitude minimum long m i n (degrees) 0.0 

Gamma-ray Intensity Skymap Longitude m cix i iiru m loTig max (degrees) 360.0 

Gamma-ray Intensity Skymap Latitude minimum lat m i n (degrees) +10.0 

Gamma-ray Intensity Skymap Latitude maximum lat max (degrees) +20.0 

Binsize in Longitude for Gamma-ray Intensity Skymaps di ong (degrees) 0.50 

Binsize in Latitude for Gamma-ray Intensity Skymaps d[ a t (degrees) 0.50 

Diffusion Coefficient Normalization Dq xx (10 28 cm 2 s _1 ) 6.1 

Diffusion Coefficient Index Below Break Rigidity D g \ 0.33 

Diffusion Coefficient Index Above Break Rigidity D g 2 0.33 

Diffusion Coefficient Break Rigidity D rigidbr (10 3 MV) 4.0 

Alfven Speed va (km s" 1 ) 30 
Nuclear Break Rigidity nuc r i g id br 

(10 3 MV) 10.0 

Nucleus Injection Index Below Break Rigidity nuc g \ 2.00 

Nucleus Injection Index Above Break Rigidity nuc g 2 2.43 

Proton Flux Normalization (10 -9 cm -2 sr" _1 s -1 MeV -1 ) 4.90 

Proton Kinetic Energy for Normalization (10 5 MeV) 1.00 

Electron Break RigidityO electron r i g idbrO (10 4 MV) 3.0 

Electron Break Rigidity electron r i g idbr (10 9 MV) 1.0 

Electron Injection Index Below Break RigidityO electron^ 2.20 

Electron Injection Index Above Break RigidityO and Below Break Rigidity electron g i 2.54 

Electron Injection Index Above Break Rigidity electron^ 2.5 

Electron Flux Normalization (10~ 10 cm" 2 sr" "is" 1 MeV" 1 ) 4.0 

Electron Kinetic Energy for Normalization (10 MeV) 3.45 
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fragmentation. 

The photon flux is proportional to the square of the number density of the DM particles 
(p/m x ) 2 , the annihilation cross section (crv) and the spectrum of photons dN/dE 1 per an- 
nihilation. The flux observed is found by integrating the number density squared along the 
line-of-sight connecting from the source to the observer, given by 

$ = V^T— 7 / dsp 2 (s^), (15) 



2 dEry Anm x j linc of sight 

where s runs along the line of sight and ip is the angle from the direction of the Galactic 
Center. Here the factor of 2 in the denominator accounts for particle-antiparticle annihilation 
since we are dealing with Dirac or complex scalar DM. We employ the isothermal profile in 



running the GALPROP [24]. 



The contribution from annihilation of DM to the observed diffuse photon spectrum is 
shown in each Fig. Q]- Fig. [3] (the long dashed curve near the bottom of each figure.) The 
position of the peak of the DM curve corresponds to about 0.1m x , as can be clearly spotted 
at each of these figures. When we add the DM contribution to the total background curve, 
we see the deviation from the data points. We can quantify the deviation from the data 
by calculating the total x 2 as a function of DM mass and the interaction strength of the 
annihilation (given by A 2 ). 

Here we adopt a simple statistical measure to quantify the effect of each operator. We 
calculate the 3a limit on each scale Aj. We assume the data agree well with the expected 
background, and then we calculate the x 2 with finite Aj's until we obtain a x 2 difference of 
Ax 2 = X 2 Xbkgd — 9 (3<t). The DM curve in each figure corresponds to a 3<r deviation 
from the data points. We calculate the limit for each operator and list them in Table IIII1 
For those unsuppressed operators the limit is of order 1 TeV for m x = 50 — 500 GeV. But for 
those operators suppressed by the velocity of the DM, light quark masses or strong coupling 
constant, the limit is significantly weaker of order 0.01 — 0.1 TeV. 

Once we obtained the lower limits on A for each DM mass, we can calculate the corre- 
sponding upper limits on (<jv)(xx ~~ QQ) as a function of DM mass. We found that these 
upper limits on (crv)(xx ~> QQ) are approximately independent of the operators with q and 
q in the final state. We understand this as the fragmentation process of q or q into ir° 
followed by ir° — > 77 decays does not depend on how the quarks are produced. As long as 
the quarks produced have the same energy, the fragmentation rate or pattern should be the 
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FERMI-LAT y-ray data ■ 

PAMELA p/p data 

_l I I I I I I I I L_ 

50 100 150 200 250 300 350 400 450 500 

m x (GeV) 

FIG. 4. The 3cr upper limits on the annihilation cross section crv(xx ~ > QQ) versus the DM 
mass due to the Fermi-LAT photon-flux data (10° < 161 < 20°, 0° < / < 360°) and due to the 

n 

PAMELA antiproton-flux data [13J. The limits are approximately independent of the operators. 
The annihilation cross sections above the curves are ruled out. The x-axis is at the value 3 x 
10 -26 cm 3 s _1 which is the required annihilation cross section to give the correct thermal relic 
density. 



same. We show in Fig. H] the upper limits on (av) ( yy — > qq) allowed by the Fermi-LAT 
photon-flux data, as well as those obtained in Ref. 13J using the antiproton-flux data of 
PAMELA . It is clear in Fig. @]that the antiproton-flux data from PAMELA does give 
a stronger constraint on the annihilation cross section (crv)(xx ~~ >* QQ) than the Fermi-LAT 
photon-flux data for lighter DM (50 — 300 GeV), while Fermi-LAT gamma-ray data does 
constrain stronger for heavier DM mass (300 — 500 GeV). 
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TABLE III. The 3er lower limits on the effective sale A of each operator listed in Table HI We take 
the coefficient C = 1 with m x = 50, 100, 200 and 500 GeV. 



Operators 




A (TeV) 






m x (GeV) 


= 50 


100 


200 


500 


Dirac DM, (axial) vector/tensor exchange 


Oi = (xYx) (qi^q) 


0.87 


1.15 


1.46 


1.94 


o 2 = (xtVx) (qi^q) 


0.025 


0.033 


0.042 


0.055 


O3 = (x7 M x) (<77m7 5 <?) 


0.87 


1.15 


1.46 


1.94 


O4 = (x7 M 7 5 x) (qi^i b q) 


0.13 


0.12 


0.11 


0.10 


o 5 = ixo^x) (q^q) 


1.04 


1.36 


1.74 


2.31 


O e = (xcr^7 5 x) (qcJ^q) 


1.04 


1.36 


1.74 


2.31 


Dirac DM, (pseudo) scalar exchange 


O7 = (xx) (qq) 


0.009 


0.011 


0.014 


0.017 


Os = (x7 5 x) (qq) 


0.094 


0.11 


0.14 


0.17 


O9 = (xx) (qj 5 q) 


0.009 


0.011 


0.014 


0.017 


O10 = (x7 5 x) (qi 5 q) 


0.094 


0.11 


0.14 


0.17 


Dirac DM, gluonic 


11 = (XX)G^ 


0.011 


0.013 


0.017 


0.024 


Oi2 = (x7 5 x) 


0.11 


0.13 


0.17 


0.24 


O13 = (XX) G^G"" 


0.011 


0.013 


0.017 


0.024 


Oi4 = (x7 5 x)G^G^ 


0.11 


0.13 


0.17 


0.24 


Complex Scalar DM, (axial) vector exchange 


O15 = (x f ^uX) (<77^) 


0.025 


0.033 


0.042 


0.055 


Oie = (X f ?x) (97^7^) 


0.025 


0.033 


0.042 


0.055 


Complex Scalar DM, (pseudo) scalar exchanj 










On = (x ] x) (qq) 


0.11 


0.10 


0.095 


0.083 


Ois = (x f x) (qi b q) 


0.11 


0.10 


0.095 


0.083 


Complex Scalar DM, gluonic 


O19 = (x ] X) G^G^ 


0.13 


0.13 


0.13 


0.14 


O 20 = (X f x) G^G^ 


0.13 


0.13 


0.13 


0.14 
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IV. DISCUSSION AND CONCLUSIONS 



Here we do a comparison with the limits obtained in collider js-5], gamma-ray lines in 
Ref. antiproton flux 13] and direct searches. 



Comparison to limits obtained using gamma-ray lines 8] against the Fermi-LAT data 
211 ] . As mentioned before, our continuum gamma-ray signals are obtained via quark 
or antiquark fragmentation, whereas in Ref. js| the discrete gamma-ray line is coming 



from one-loop dressing of t 
3a lower limit while Ref 



re effective operators. Note also that we are calculating the 
reported the 95% C.L. lower limits, which are approxi- 
mately 2a, so our limits are slightly more conservative. We found that the limits for 
operator Oi 3 are substantially better than the corresponding operators D 5 j of Ref. 
Our limits are 0.9 — 1.9 TeV for m x = 50—500 GeV, while their limits are 0.12—0.6 TeV 
for m x = 50 — 200 GeV. We see in this case (Dirac DM with vector-boson exchange), 
the contribution to photon flux through the continuum (q — > 7T° — > 7) is substantially 
better than the discrete line spectrum (via the loop process) when compared against 
the Fermi-LAT data. The limits for operators Oj^iq (Dirac DM with scalar-boson 
exchange) and On^% (complex scalar DM with scalar-boson exchange) are about the 
same as -Di_4 and Cx s 2, respectively, of Ref. 8[. On the other hand, limits for opera- 
tors Oi5,i6 (complex scalar DM with vector-boson exchange) are weaker than 63,4 of 
Ref. [8]. 



Comparison to limits obtained using antiproton flux 



ton data 25]. The approach that we did in Ref. 13] is very similar to the work here, 



13] against the PAMELA antipro- 



except that we require to see photon instead of antiproton in the final state and we are 
comparing with two entirely different categories of data. We found that at the lower 
end m x ~ 50 — 300 GeV the limits from antiproton data are somewhat better than 
those from photon-flux data. This is because of the nature of the fragmentation and 
decay chain: one fragments into antiproton directly while the other one fragments into 
7T° then followed by decays of 7r° — > 77. The peak of the antiproton spectrum occurs 
at larger energies than the peak of the photon spectrum for the same DM mass. On 
the other hand, when m x increases to about 300 — 500 GeV, limits from gamma-ray 
flux constrains stronger on the DM interactions. 
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Comparison to limits obtained using collider data 4|, |5|. The collider limits on A come 
from monojet plus missing energy data. Basically, for operators that are velocity 
suppressed, e.g., 2 ,4 among Oi_ 4 , the limits from photon flux are not as good as 
those from colliders; and vice versa for operators that are not velocity suppressed, 
e.g., Oi ; 3 among Oi_ 4 , The collider limits for Oi_4 are 200 — 300 GeV, while the limits 
from photon flux are 0.9 — 1.9 TeV for but only 25 — 100 GeV for 02,4- The 
operators for which photon flux gives better limits are Oi,3, 0$^, Og,io, Oi2,i4> Oi7,i8 
and Oi9,2o- 



Comparison to limits obtained by the direct detection experiments [26|, |27jj. Two most 
stringent experiments on the spin-independent (SI) cross sections come from CDMS 



26| and XENON100 |27j. The best upper limit from CDMS |26| is a 



10 



-li 



cm at m 



x 



70 GeV and from XENON100 



SI 

x^ 



3.8 x 



221 ^ 4n 



0.7 x 10 



-ii 



at 



50 GeV. Note that only the operators 0\ and Or for Dirac DM contribute to 
spin- independent cross sections. For 7 = m 9 (xx) (<?<?) /A 3 the SI cross section is given 
by 



71 



G 



N\2 



(16) 



where /i x tv is the reduced mass of the DM and nucleon, and 



\ A 3 



(N\qq\N) = — x 



fN 
JTq 

2 j?N 



light quarks 



27 fjlg '■ heavy quarks 



(17) 



Using the default values of f^ q and fjf adopted in DarkSUSY [28| and averaging 

si 

r x iV> 



between the neutron and proton for <Ty N , we obtain 



o* - ^(0.3769) 2 . 



(18) 



Using the limit a^ N < 10~ 44 cm 2 for m x ~ 50 - 200 GeV [27J, we obtain 



A 7 > 330 GeV . 



(19) 



This limit is much better than the limit on Oj shown in Table IIHI On the other hand, 
with the operator 0\ = (x7 M x) (97^) / ' ^\ we obtain 



a 



SI _ I^XN 

~ 256vr 
18 



h 12 



>N\ 



(20) 



where ~ + a d)i which is obtained by taking the average of neutrons and 

protons inside a nucleon, and a^ d = l/A 2 are the coefficients in front of the operators 
for the valence u and d quark. We therefore obtain 

a * N ~ 2567T Af ' W 
which gives the limit on A 1; with a^ N < 10~ 44 cm 2 , as 

Ai > 4.4 TeV . (22) 

It is about 2 — 3 times better than the corresponding limits on 0\ shown in Table IIII1 

In summary, we have used an effective interaction approach to investigate the effects 
of dark matter interactions with light quarks on diffuse photon flux coming from the mid- 
latitude (10° < |6| < 20°, 0° < / < 360°) region. We have assumed a standard halo 
density while using the GALPROP to calculate the resulting diffuse gamma-ray spectrum. 
The diffuse gamma-ray background includes tt° decays produced by scattering of cosmic 
rays with the interstellar medium, inverse Compton scattering, bremsstrahlung, synchrotron 
radiation and extragalactic diffuse gamma rays. The background can reproduce the data 
recorded by Fermi-LAT very well. Dark matter annihilation provides yet another source of 
diffuse gamma rays. The dominant mechanism is annihilation into light quarks, followed by 
fragmentation into neutral pions, which then further decay into photons. The gamma-ray 
spectrum has no particular feature but a continuum; nevertheless, the flux is large. We have 
shown that the effective interactions of the DM can give rise to a nontrivial contribution to 
the diffuse gamma-ray spectrum, as presented in various figures. 

We have successfully used the data to obtain 3a limits on the scale Aj. The best limits 
are from the Dirac DM with vector or tensor boson exchanges. The limits for and 

6 are about 1 — 2 TeV, while other operators suppressed either by the velocity of dark 
matter, light quark masses or strong coupling constant give milder limits. Note that these 
limits from photon flux are lower limits on Aj. We found that the limits obtained are very 



comparable to t 



rose obtained using collider {4,5], gamma-ray line search js] and anti-matter 



experiments [13J. However, the two operators (xx)(qq) an d (XT^xX^T/^)) which contribute 
to spin-independent cross sections, are constrained more severely by the recent XENON100 
data. 



19 



ACKNOWLEDGMENTS 



This work was supported in parts by the National Science Council of Taiwan under Grant 
Nos. 99-21 12-M-007-005-MY3 and 98-2112-M-001-014-MY3 as well as the WCU program 
through the KOSEF funded by the MEST (R31-2008-000-10057-0). 



[1] E. Komatsu et al. [WMAP Collaboration], "Seven- Year Wilkinson Microwave Anisotropy 
Probe (WMAP) Observations: Cosmological Interpretation," Astrophys. J. Suppl. 192, 18 
(2011) [arXiv: 1001 .45381 [astro-ph.CO]]. 

[2] G. Bertone, D. Hooper and J. Silk, "Particle dark matter: Evidence, candidates and con- 
straints," Phys. Rept. 405, 279 (2005) [arXiv:hep-ph/0404175| . 

[3] Q. H. Cao, C. R. Chen, C. S. Li and H. Zhang, "Effective Dark Matter Model: Relic density, 
CDMS II, Fermi LAT and LHC," i arXiv:0912.45lTl [hep-ph]. 

[4] Y. Bai, P. J. Fox and R. Harnik, "The Tevatron at the Frontier of Dark Matter Direct 
Detection," larXiv: 1005.37971 [hep-ph]. 

[5] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. P. Tait and H. -B. Yu, "Constraints 
on Dark Matter from Colliders," Phys. Rev. D82, 116010 (2010) [ arXiv: 1008. 17831 [hep-ph]]. 

[6] M. Beltran, D. Hooper, E. W. Kolb, Z. A. C. Krusberg and T. M. P. Tait, "Maverick dark 
matter at colliders," JHEP 1009, 037 (2010) |arXiv: 1002.41371 [hep-ph]]. 

[7] G. D. Mack, T. D. Jacques, J. F. Beacom, N. F. Bell and H. Yuksel, "Conservative Con- 
straints on Dark Matter Annihilation into Gamma Rays," Phys. Rev. D 78, 063542 (2008) 
|arXiv:0803.0T57l [astro-ph]] . 

[8] J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. P. Tait and H. B. P. Yu, "Gamma 
Ray Line Constraints on Effective Theories of Dark Matter," larXiv:10 09.0008 [hep-ph]. 

[9] D. T. Cumberbatch, Y. L. Tsai and L. Roszkowski, "The impact of propagation uncertainties 
on the potential Dark Matter contribution to the Fermi LAT mid-latitude gamma-ray data," 
Phys. Rev. D 82, 103521 (2010) |arXiv: 1003 .28081 [astro-ph.HE]]. 
[10] D. Hooper and L. Goodenough, "Dark Matter Annihilation in The Galactic Center As Seen 
by the Fermi Gamma Ray Space Telescope," Phys. Lett. B 697, 412 (2011) [ arXiv: 1010.27521 
[hep-ph]]. 

20 



[11] J. Fan, M. Reece and L. T. Wang, "Non-relativistic effective theory of dark matter direct 

detection," larXiv:1008.159Tl [hep-ph]. 
[12] K. Cheung, K. Mawatari, E. Senaha, P. Y. Tseng and T. C. Yuan, "The Top Window for 

Dark Matter," JHEP 1010, 081 (2010) [arXiv:1009.06l8l [hep-ph]]. 
[13] K. Cheung, P. Y. Tseng and T. C. Yuan, "Cosmic Antiproton Constraints on Effective Inter- 
actions of the Dark Matter," JCAP 1101, 004 (2011) [arXiv:1011.23l0l [hep-ph]]. 
[14] C. Boehm, T. Delahaye, P. Salati, F. Staub and R. K. Singh, "Implication of the PAMELA 

antiproton data for dark matter indirect detection at LHC," JCAP 1006, 013 (2010) 

|arXiv:0907.45lTl [hep-ph]]. 
[15] K. Cheung, J. Song and P. Y. Tseng, "Cosmic positron and antiproton constraints on the 

gauge-Higgs Dark Matter," JCAP 1009, 023 (2010) |arXiv:1007.0282l [hep-ph]]. 
[16] P. Scott, J. Conrad, J. Edsjo, L. Bergstrom, C. Farnier and Y. Akrami, "Direct Constraints 

on Minimal Supersymmetry from Fermi-LAT Observations of the Dwarf Galaxy Segue 1," 

JCAP 1001, 031 (2010) |arXiv:0909.3300l [astro-ph.CO]]. 
[17] D. Hooper and L. Goodenough, "Dark Matter Annihilation in The Galactic Center As Seen 

by the Fermi Gamma Ray Space Telescope," Phys. Lett. B 697, 412 (2011) [ arXiv: 1010.27521 

[hep-ph]]. 

[18] P. Ciafaloni, M. Cirelli, D. Comelli, A. De Simone, A. Riotto and A. Urbano, "On 
the Importance of Electroweak Corrections for Majorana Dark Matter Indirect Detection," 
larXiv: 1104^9961 [hep-ph]; P. Ciafaloni, D. Comelli, A. Riotto, F. Sala, A. Strumia and A. Ur- 
bano, "Weak Corrections are Relevant for Dark Matter Indirect Detection," JCAP 1103, 019 
(2011) [arXiv: 1009.02241 [hep-ph]]. 

[19] N. F. Bell, J. B. Dent, A. J. Galea, T. D. Jacques, L. M. Krauss and T. J. Weiler 
"W/Z Bremsstrahlung as the Dominant Annihilation Channel for Dark Matter, Revisited," 
larXiv: 1104.38231 [hep-ph]. 

[20] K. Nakamura et al. (Particle Data Group), J. Phys. G37, 075021 (2010). 

[21] A. A. Abdo et al. [The Fermi-LAT Collaboration], "The Spectrum of the Isotropic Diffuse 
Gamma-Ray Emission Derived From First- Year Fermi Large Area Telescope Data," Phys. 
Rev. Lett. 104, 101101 (2010) [arXiv: 1002.36031 [astro-ph.HE]]. 

[22] T. Sjostrand, S. Mrenna and P. Z. Skands, "A Brief Introduction to PYTHIA 8.1," Comput. 
Phys. Commun. 178, 852-867 (2008) [arXiv:0710.3820l [hep-ph]]. 



21 



[23] S. Albino, B. A. Kniehl and G. Kramer, "Fragmentation functions for light charged 
hadrons with complete quark flavour separation," Nucl. Phys. B 725, 181 (2005) 
|arXiv:hep-ph/0502188] . 

[24] A. W. Strong, I. V. Moskalenko, T. A. Porter, G. Johannesson, E. Orlando and S. W. Digel, 

"The GALPROP Cosmic-Ray Propagation Code." la?Xlv:0907.0559l [astro-ph.HE]. 
[25] O. Adriani et al. [PAMELA Collaboration], "PAMELA results on the cosmic-ray antipro- 

ton flux from 60 MeV to 180 GeV in kinetic energy," Phys. Rev. Lett. 105, 121101 (2010) 

|arXiv:1007.082T1 [astro-ph.HE]]. 
[26] Z. Ahmed et al. [The CDMS-II Collaboration], Science 327, 1619-1621 (2010). 

|arXiv:0912.3592l [astro-ph.CO]]. 
[27] E. Aprile et al. [XENON100 Collaboration] , |ar2Gyin042549J [astro-ph.CO]. 
[28] P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schelke and E. A. Baltz, JCAP 0407, 008 

(2004) |arXiv: astro-ph/0406204| . 



22 



